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INTRODUCnOR 

SHADOW  Is  a  l6o4  COP-controUed  program  designed  to  confute  the  times  at 
which  a  satellite  enters  and  leaves  the  shadow  of  the  earth. 

1.0  IDEHTmCATION 

Congpute  shadow  times  (SHADOW) 

2.0  IVNCnON 

Ibe  Inputs  to  this  program  are  a  reset  tape  containing  orblteil  conditions 
and  a  function  card  specifying  a  time  range  and  a  time  Increment.  Obe  latter 
Is  the  time  between  cosgputatlons  and  Is  equivalent  to  accuracy  desired.  The 
outinit  Is  an  on-line  printout  of  shadow  periods  (see  Section  6.0). 

3.0  MATHEMATICAL  EOBMUIATION 

The  following  assungptlons  are  made: 

a.  The  shadow  of  the  earth  Is  cylindrical  and  non-fuzzy. 

b.  The  positions  of  the  satellite  and  the  sun  are  known  for  a  given  time, 
t.  In  lnez*tlal  coordinates. 

The  coordinates  of  the  sun  at  time,  t,  are  (x,y,z).  The  axis  of  the 
cylinder  formed  by  the  shadow  of  the  earth  passes  through  the  point  (-x,-y,-z). 

The  coordinates  of  the  satellite  at  time,  t,  are  ^  )* 

IX  I  -Z 

The  equation^ defines  an  earth-centered  sphere  of  radius, |0  , 
with  the  satellite  on  the  svirface. 
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The  axis  of  the  cylindrical  shadow  pierces  this  sphere  at  x',  y',  z', 


i*ere 


X'  -  -  X 

y'  • 

z'  Y 


Graphically^  the  situation  can  be  represented  as  follows : 
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l!be  distance,  d,  from  (x’,  y',  z' )  to  the  clrcxilar  Intersection  of  the 
cycUnder  and  the  sphere  Is  given  by: 

d^  ■  ~  where  r  ■  earth  radius 

The  distance,  D>  from  (x',  y',  z')  to  the  satellite  is  given  by 
T?  -  (x'  -S^f  +  (y  -rif  +  (z* 

If  0^>d^,  the  satellite  is  in  sunlight. 

2  2 

If  D  4  d  ,  the  satellite  is  in  shadow. 

Please  note  that  the  analysis  above  does  not  assume  a  circular  or  near- 
circular  orbit .  A  new  p  is  conputed  for  each  t . 

4.0  CALLING  SEQUENCE 

This  program  is  normally  called  by  a  function  card  with  the  format: 


*  SHADOW  A  B  C  D  E  F.F  0  H  I  J  K.K  L.L 

A  Vehicle  number 
B  Month 
C  Day 

D  Hotir  Stcurt  Time 

E  Minute 
F.F  Second 
G  Month 
H  Day 

I  Hour  Stop  Time 

J  Minute 

K. K  Second 

L. L  Confutation  interval  in  seconds 

An  option  exists  to  specify  start  and  stop  times  by  rev  nuniber.  The 
function  card  has  a  format: 
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♦  SHADOW  A  B.  C.  D.D 


A  Vehicle  Nvodber 

B  Starting  Rev  (decimal  point  mandatory) 

C  Qidlng  Rev  (decimal  point  mandatory) 

D.D  Computation  Interval  In  seconds 

A  fractional  rev  number  may  be  Input  vlthout  harming  the  program,  but 
the  fractional  part  will  be  ignored.  For  consistency  between  the  function 
card  and  printed  output,  It  is  recommended  that  zero,  or  blank,  follow  the 
decimal  point  In  the  rev  nvunber. 

Confutation  Is  carried  out  to,  but  not  Including,  the  terminal  rev. 

Ihe  last  field  on  the  card,  computation  Interval,  controls  the  accuracy 
of  the  printout  and  the  time  consumed  by  the  program.  For  example.  If  the 
cosputatlon  interval,  were  60,  the  printout  would  be  accurate  to  within  one 
minute.  If  the  cosputatlon  Interval  were  1,  the  printout  accuracy  would 
be  one  second  and  the  program  would  take  roughly  sixty  times  longer.  Crude 
timing  estimates  Indicate  that  for  a  computation  Interval  of  sixty  seconds, 
the  program  requires  between  one  and  two  seconds  per  rev. 

^.0  DEmiLED  PROCESSES 

The  program  conputes  the  inertial  coordinates  of  the  satellite  at  time, 
t,  by  successive  branches  to  the  subroutines  TTS  and  FTR.  nie  coordinates 
of  the  sun,  obtained  by  table  look-up,  are  taken  from  the  American  Ephemerls 
and  Hautlce^.  Almanac  and  contain  one  value  of  x,  y,  and  z  for  each  day  of 
the  year  in  the  following  format. 

Location  Contents 

L  X 

L+1  y 

h*2 


z 


Deceniber  31>  preceding  yesur 
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Location 


Contents 


L+3 

X 

y 

Jamuary  1,  current  year 

L+5 

z 

L+1095 

* 

X 

L+1096 

y 

December  31 »  current  yearn 

I/<-l097 

z 

L+109e 

X 

L+1099 

y 

Jamuaxy  1,  next  year 

L+noo 

z 

^  get  the  sun  coordinates  at  a  given  time 

t,  the  program  performs  a 

linear  Interpolation  hetveen  successive  triads 

of  tabuleur  (x,  y,  z). 

Ihe  entries  In  the  table 

cam  be  tramscrlbed  directly  from  the  almanac. 

Assuming  the  first  address  of  the  table  to  be  tagged  SUNOOORD,  the  first  few 

instruction  cards  for  the  1963  version  of 

the  table  would  be: 

1 

2 

4 

Columns  123U56789032 

OI23U5678 

1  (Comments  field) 

SUNCOORD  DEC 

+.15096i«8 

X,  Dec  31 

DEC 

-.8911*589 

y,  Dec  31 

DEC 

-.38658IH 

z,  Dec  31 

DEC 

+.l6824Jf5 

X,  Jam  1 

DEC 

-.8888335 

y,  Jan  1 

DEC 

-.3951*1*82 

z,  Jan  1 

(etc. ) 
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6.0  OUTPUTS 

Every  time  the  satellite  leaves  the  shadow  of  the  eso^th,  the  program 
prints  on  line  the  date,  time,  rev  niiniber,  latitude  (geodetic),  and  longitude 
at  \dilch  It  enterred  and  left.  When  the  stop  time  is  reached,  the  program 
prints  out  the  percentage  of  time  In  shadow  (l.e.,  the  sum  of  all  the  shaded 
periods  divided  hy  the  time  range  on  the  function  card).  A  sample  printout  Is 
shown  In  Appendix  A. 

7.0  IMTERFACES  AND  RESTRICTIONS 

Ihe  program  usee  the  following  subroutines. 

Code  name 

RESET 

PTR 

LTT 

TTE 

SETUP 

SUBERR 

RADDEG 

OUTPUT 

FLOAT 

FIX 

RALONG 

UTMAC 

MACGUT 

SORT 

ATANl 

TAN 

Once  a  year  the  table  of  sun  coordinates  (see  Section  5«0)  must  be  updated, 
and  the  program  reassembled. 

Reset  tapes  containing  orbital  conditions  must  be  mounted  on  tape  drives 
2,  9>  and  10. 


t 
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To  start  from  scratch,  the  follovdng  function  sequence  Is  reconmended : 
*IRT  ... 


*WNRr  ... 

*INJ]fTJN... 


(Each  card  conq>leted  appropriately) 


♦SHADOW  ... 


f 
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8.0  FLOW  CHAKT 


£3iTER 


SAVE  Index  Registers,  Print  Headers 
Read  File  1  Reset  ^pe.  Run  SEIUP 


Set  Start  cmd  Stop 


Bines  from  function  card 


COMPUTE  Compute  eluents  this  t,  polar  and  rectangular 
Cong^nite  day  of  the  year 
Compute  ^niATlO  ■  Seconds  Into  this  day/86U00 
Interpolate  to  get  sun  coordinates 
Compute  Ra^e  of  Sun 


Con^mte  x'y  y',  z*,  point  where  axis  of  cylindrical  shadow 
Pierces  s^ere  containing  satellite 

Compute  d  ,  distance  from  x',  y',  z'  to  circumference  of 
shaded  cap 


A 


Con^te  Vr,  distance  from  x',  y',  z’  to  satellite 

1)2  _^2  £  0^  _Ye6 - ^  SHADE  This  means  we  are  In  shadow 

I  Were  we  In  shadow  last  time? 

KO  (SHFIA6  «  1)7 


r 


This  means  satellite  In 
sunshine.  Were  we  In 

sunshine  last  time?  -YES  - 

(SHFIAG  -  0?) 

I 

NO 

>1 

This  means  we  Just  entered 
svinshlne.  Print  out  one  line 
of  data 

4. 

Set  SHFLAO  -  0 

i 

Accumulate  total  shadow 
time, 


YES 

4 


NO 

i 

Save  entering  times 
Set  SHPIAO  -  1 


SHP^C 


NElCr  Increment  TAU 

I 

Have  we  exceeded  stop  time 
(TAU-WHISTLB  >  O)? 


to 


O^PUTE 


YES 

Reset  Index  Register 
Print  Summau^y 
Return  to  OOP 


APPEHDIX 
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